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Abstract 

The Monte Carlo simulation of Majorana fermions is discussed 
on the example of supersymmetric Yang-Mills (SYM) theory. 



1 Introduction 

Majorana fermions play an important role in supersymmetric quantum 
field theories which are the favorite candidates for the extension of the 
Standard Model of elementary particle interactions beyond the presently 
available energy range. It is generally assumed that the scale where su- 
persymmetry becomes manifest is near to the presently explored elec- 
troweak scale and that the supersymmetry breaking is spontaneous. An 
attractive possibility for spontaneous supersymmetry breaking is to ex- 
ploit non-perturbative mechanisms in supersymmetric gauge theories. 

Non-perturbative features of supersymmetric theories can be derived 
analytically, as for instance the basic work of Seiberg and Witten |L] shows, 
and can also be obtained by numerical Monte Carlo simulations on a 
lattice. The numerical approach requires the introduction of Majorana 
fermion fields on a four- dimensional Euclidean lattice in space and imag- 
inary time. 

The simplest supersymmetric gauge theory is the supersymmetric ex- 
tension of Yang-Mills gauge theory. It is the gauge theory of a massless 
Majorana fermion, called "gaugino", in the adjoint representation of the 
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gauge group. The Euclidean action density of a gauge theory in the adjoint 
representation can be written as 

£ = \F^F^ + ^T^(p^) a + m s r\ a . (1) 

Here F£ v denotes the field strength tensor and A a is the Grassmannian 
fermion field, both with the adjoint representation index a. rrig is the 
gaugino mass which has to be set equal to zero for supersymmetry. For a 
Majorana fermion A a and A are not independent but satisfy 

A = X T C , (2) 

with C the charge conjugation Dirac matrix. This definition is based 
on the analytic continuation of Green's functions from Minkowski to Eu- 
clidean space 0. 

2 Lattice formulation 

In order to define the path integral for a Yang Mills theory with Majorana 
fermions in the adjoint representation, let us first consider the familiar 
case of Dirac fermions 0. (For a general reference on lattice quantum 
field theory see this book.) Let us denote the Grassmanian fermion fields 
in the adjoint representation by ^ and ip a x . Here Dirac spinor indices are 
omitted for simplicity and a stands for the adjoint representation index 
(a = 1, .., N% — 1 for SU(iV c ) ). The fermionic part of the Wilson lattice 
action is 

s f = E{C^ 

X 

- K £ [C+a Wl + + - 7m)^. +A ] } • ( 3 ) 

Here K is the hopping parameter, the Wilson parameter removing the 
fermion doublers in the continuum limit is fixed to r = 1 and the matrix 
for the gauge-field link in the adjoint representation V Xfl is defined from 
the fundamental link variables U Xfl according to 

= V ab ^[U] = 2Tr(UlT a U x ,T b ) = V; b ^ = V~ h ^ . (4) 
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The generators T a = ^X a satisfy the usual normalization Tr (A A&) = \. 
In the simplest case of SU(2) (N c = 2) we have, of course, T a = \r a with 
the isospin Pauli- matrices r a . The normalization of the fermion fields in 
(|3p is the usual one for numerical simulations. The full lattice action is 
the sum of the pure gauge part and fermionic part: 

S = S g + S f . (5) 

The standard Wilson action for the SU(iV c ) gauge field S g is a sum over 
the plaquettes 

pi v iv c J 

with the bare gauge coupling given by (3 = 2N c /g 2 . 

In order to obtain the lattice formulation of a theory with Majorana 
fermions let us note that out of a Dirac fermion field it is possible to 
construct two Majorana fields: 

A« = ity + C^p T ) , A< 2 > = + Ct) (7) 

with the charge conjugation matrix C. These satisfy the Majorana con- 
dition 

\ U) = \^ T C (j = 1,2). (8) 

The inverse relation of (|7|) is 

V> = ^(A (1) + *A( 2 )) , ^ = cf = -L(A« - • (9) 

In terms of the two Majorana fields the fermion action Sf in eq. (Q) can 
be written as 

2 



Z x- j=l 

- k j: te; wi + + A?v a ^(i - 7 ,)Aifj }(io) 



^=1 

For later purposes it is convenient to introduce the fermion matrix 

Qyd,xc = Qyd,xc\U~\ = ^yx^dc 
4 

- [s y ,x +fi (i + ^)v dc ^ + w(!-7Xj • ( n ) 
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Here, as usual, fi denotes the unit vector in direction \x. In terms of Q we 
have 

s f = E tyQv^i = ^EE Xf'^/.i.Ar , (12) 

xc,yd j=l xc,yd 

and the fermionic path integral can be written as 

2 

/ [#rfV]e" S/ = / [##]e^ w = detQ = n / [dA^e"^ . (13) 

This shows that the path integral over the Dirac fermion is the square of 
the path integral over the Majorana fermion and therefore 

J [d\]e-^ QX = i^detQ . (14) 

As one can see here, for Majorana fields the path integral involves only 
[dA^] because of the Majorana condition in @. 

The relation (|14]) leaves the sign on the righ hand side undetermined. 
A unique definition of the path integral over a Majoran fermion field, 
including the sign, is given by 

J[d\]e-^ QX = J[dX}e^ XMX = Pf(M) (15) 

where M is the antisymmetric matrix defined as 

M = CQ = -M T . (16) 

The square root of the determinant in eq. fll4"D is a Pfaffian H. This 
can be defined for a general complex antisymmetric matrix M a p = —Mp a 
with an even number of dimensions (1 < a, /3 < 2N) by a Grassmann 
integral as 

Pf (M) = / [d4]e-fa">*+» = j^e aih ... aNPN M aiPl . . . M QjvAv . (17) 

Here, of course, [d<p] = d(p 2 N ■ ■ -d&ii an d e is the totally antisymmetric 
unit tensor. 

It is now clear that the fermion action for a Majorana fermion in the 
adjoint representation A" can be defined by 
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This together with (|5])-(|6]) gives a lattice action for the gauge theory of 
Majorana fermion in the adjoint representation. In order to achieve super- 
symmetry one has to tune the hopping parameter (bare mass parameter) 
K to the critical value K cr (f3) in such a way that the mass of the fermion 
becomes zero. 

The path integral over A is defined by the Pfaflian Pf(CQ) = Pf(M). 
By this definition the sign on the right hand side of eq. QUI) is uniquely 
determined. The determinant det(Q) is real because the fermion matrix 
in (|ll|) satisfies 

Q f = 75Q75 , Q = 75Q = Q ] • (19) 

Moreover one can prove that det(Q) = det(Q) is always non-negative. 
This follows from the relations 

CQC- 1 = Q T , BQB- 1 = Q T , (20) 

with the charge conjugation matrix C and B = C75. It follows that every 
eigenvalue of Q and Q is (at least) doubly degenerate. Therefore, with 
the real eigenvalues Aj of the Hermitean fermion matrix Q, we have 

det(Q) = det(Q) = JJ \? > 0. (21) 

i 

Since according to the above discussion 

det(Q) = det(M) = [Pf(M)] 2 , (22) 
the Pfaflian Pf(M) has to be real - but it can have any sign. 



3 Numerical simulations of SYM theories 



In order to perform Monte Carlo simulations of SYM theory one needs a 
positive measure on the gauge field which allows for importance sampling 
of the path integral. Therefore the sign of the Pfaflian can only be taken 



into account by reweighting. According to (22) the absolute value of the 
Pfaffian is the non-negative square root of the determinant therefore the 
effective gauge field action is 0: 

Scv=PY,( 1 - 7p U p) - \ lo § det QM ■ ( 23 ) 
pi ~ / _ 
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The factor | in front of log det Q shows that we effectively have a flavour 
number Nf — | of adjoint fermions. The omitted sign of the Pfaffian can 
be taken into account by reweighting the expectation values according to 

_ (AsignPf(M)) cy 
{A) ~ (signPf(M)) cy ' ^ 

where (. . .)cv denotes expectation values with respect to the effective 
gauge action Scv- This may give rise to a sign problem which will be 



discussed in section 3.1 



The fractional power of the determinant corresponding to (p3|) can be 
reproduced, for instance, by the hybrid molecular dynamics algorithm Jj] 
which is, however, a finite step size algorithm where the step size has to be 
extrapolated to zero. An "exact" algorithm where the step size extrapola- 
tion is absent is the two-step multi-bo sonic (TSMB) algorithm |7], ||. The 
first large scale numerical simulation of SYM theory has recently been 
performed by the DESY-Miinster-Roma collaboration using the TSMB 
algorithm f| llOL HI|. 



3.1 The "sign problem" 

The Pfaffian resulting from the Grassmannian path integrals for Majorana 
fermions (|T5|) is an object similar to a determinant but less often used. 
As shown by (|T7D, Pf(M) is a polynomial of the matrix elements of the 
2iV-dimensional antisymmetric matrix M = —M T . Basic relations are [f| 

M = P T JP, Pf (M) = det(P) , (25) 

where J is a block-diagonal matrix containing on the diagonal 2® 2 blocks 
equal to e = io<i and otherwise zeros. Let us note that from these relations 
the second equality in eq. ( P2| ) immediately follows. 

The form of M required in fl25|) can be achieved by a procedure anal- 
ogous to the Gram-Schmidt orthogonalization and, by construction, P is 
a triangular matrix (see [[HJ). This gives a numerical procedure for the 
computation of P and the determinant of P gives, according to (^5|), the 
Pfaffian Pf(M). Since P is triangular, the calculation of det(P) is, of 
course, trivial. 

This procedure can be used for a numerical determination of the Pfaf- 
fian on small lattices. On lattices larger than, say, 4 3 ■ 8 the computation 
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becomes cumbersome due to the large storage requirements. This is be- 
cause one has to store a full Q ® Q matrix, with Q being the number 
of lattice points multiplied by the number of spinor-colour indices (equal 
to 4(iV^ — 1) for the adjoint representation of SU(iV c )). The difficulty of 
computation is similar to a computation of the determinant of Q with 
LU -decomposition. 

Fortunately, in order to obtain the sign of the Pfaffian occurring in 
the reweighting formula (|23|), one can proceed without a full calculation 
of the value of the Pfaffian. The method is to monitor the sign changes of 
Pf(M) as a function of the hopping parameter K. According to fl2"T|), the 
hermitean fermion matrix for the gaugino Q has doubly degenerate real 
eigenvalues therefore 

Q/2 

det M = det Q = f[ A? , (26) 
i=i 

where Aj denotes the eigenvalues of Q. This implies 

a/2 n/2 _ 

|Pf(M)| = niAi| , =^ Pf(M) = nAV (27) 
i=l 1=1 



The first equality trivially follows from (|22|). The second one is the con- 
sequence of the fact that Pf (M) is a polynomial in K which cannot have 
discontinuities in any of its derivatives. Therefore if, as a function of K, 
an eigenvalue Aj (or any odd number of them) changes sign the sign of 
Pf(M) has to change, too. Since at K = we have Pf(M) = 1, the num- 
ber of sign changes between K = and the actual value of K, where the 
dynamical fermion simulation is performed, determines the sign of Pf (M). 
This means that one has to determined the flow of the eigenvalues of Q 
through zero [|T2]]. 

The spectral flow method is well suited for the calculation of the sign 
of the Pfaffian in SYM theory. An important question is the frequency 
and the effects of configurations with negative sign. A strongly fluctuating 
Pfaffian sign is a potential danger for the effectiveness of the Monte Carlo 
simulation because cancellations can occur resulting in an unacceptable 
increase of statistical errors. The experience of the DESY-Miinster Collab- 
oration shows, however, that below the critical line K cr ((3) corresponding 
to zero gaugino mass (trig = 0) negative Pfaffians practically never appear 
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[jrO| , pTl|| . Above the critical line several configurations with negative Pfaf- 
fian have been observed but their role has not yet been cleared up to now. 
Since supersymmetry is expected to be realized in the continuum limit at 
wig = 0, the negative signs of the Pfaffian can be avoided if one takes the 
zero gaugino mass limit from nig > corresponding to K < K cr . In this 
sense there is no "sign problem" in SYM which would prevente a Monte 
Carlo investigation. 

The presence or absence of negative Pfaffians in a sample of gauge 
configurations produced in Monte Carlo simulations can be easily seen 
even without the application of the spectral flow method. In case of sign 
changes the distribution of the smallest eigenvalues of the squared fermion 
matrix Q 2 shows a pronounced tail reaching down to zero |L3| . The ab- 
sence of a tail shows that there are no negative Pfaffians. 

Concerning this "sign problem" let us note that a very similar phe- 
nomenon appears also in QCD because the Wilson-Dirac determinant of a 
single quark flavour can also have a negative sign. Under certain circum- 
stances the sign of the quark determinant plays an important role. This 
is the case, for instance, at large quark chemical potential in a QCD-like 
model with SU(2) colour and staggered quarks in the adjoint representa- 
tion which has recently been studied by the DESY-Swansea Collaboration 
HJ. This investigation also revealed an interesting feature of the TSMB 



algorithm, namely its ability to easily change the sign of eigenvalues of the 
hermitean fermion matrix (and hence the sign of the determinant or Pfaf- 
fian). This is in contrast to algorithms based on finite difference molecular 
dynamics equations as, for instance, the HMD || algorithms. 



4 Outlook 

Numerical simulations of N — 1 supersymmetric theories require to deal 
with Majorana fermions on a Euclidean lattice. Such simulations are 
feasible with presently available computer technology and using well es- 
tablished simulation algorithms - at least in the relatively simple case of 
supersymmetric Yang-Mills theories || |10, [11], [14], [lf|]. Many interest- 



ing questions are waiting for detailed answers. Just to mention a few: 
the behaviour of the phase transition at the supersymmetric point, the 
spectroscopy of supersymmetric multiplets in the particle spectrum and 
the Ward-Takahashi identities proving the realization of supersymmetry 
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in quantum field theories. 
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